arXiv:1506.05459v2 [astro-ph.GA] 29 Jul 2015 


Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 30 July 2015 (MN KTeX style file v2.2) 


ALMA and Herschel reveal that X-ray selected AGN and 
main-sequence galaxies have different star formation rate 
distributions 


J. R. Mullaney^*, D. M. Alexander^, J. Aird^, E. Bemhard\ E. Daddi"^, A. Del Moro^, 
M. Dickinson^, D. Elbaz"^, C. M. Harrison^, S. Juneau^, D. Liu^, M. Pannella^, 

D. Rosario^, R Santini^, M. Sargent^, C. Schreiber"^, J. Simpson^, E. Stanley^ 

^ Department of Physics and Astronomy, The University of Sheffield, Hounsfield Road, Sheffield, S3 7RH, UK 
^Centre of Extragalactic Astronomy, Department of Physics, Durham University, South Road, Durham, DHl 3LE, UK 
^Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CBS OHA, UK 

‘^Laboratoire AIM, CEA/DSM-CNRS-Universite Paris Diderot, Irfu/Service dAstrophysique, CEA-Saclay, Orme des Merisiers, 91191 Gif-sur-Yvette, France 
^National Optical Astronomy Observatories, 950 N Cherry Avenue, Tucson, AZ 85719, USA 
^ Universitdts-Sternwarte Miinchen, Scheinerstr. 1, D-81679 Miinchen 

Max-Planck-Institut fr Extraterrestrische Physik (MPE), Postfach 1312, 85741, Garching, Germany 
^INAF-Osservatorio Astronomico di Roma, via di Frascati 33,1-00040 Monte Porzio Catone, Roma, Italy 
'^Astronomy Centre, Department of Physics and Astronomy, University of Sussex, Brighton, BNl 9Q11, UK 


Date Accepted 


ABSTRACT 

Using deep Herschel and ALMA observations, we investigate the star formation rate 
(SFR) distributions of X-ray selected AGN host galaxies at 0.5 < z < 1.5 and 1.5 < z < 4, 
comparing them to that of normal, star-forming (i.e., “main-sequence”, or MS) galaxies. We 
find 34-55 per cent of AGNs in our sample have SFRs at least a factor of two below that of the 
average MS galaxy, compared to ft! 15 per cent of all MS galaxies, suggesting significantly 
different SFR distributions. Indeed, when both are modelled as log-normal distributions, the 
mass and redshift-normalised SFR distributions of X-ray AGNs are roughly twice as broad, 
and peak ft 0.4 dex lower, than that of MS galaxies. However, like MS galaxies, the nor¬ 
malised SFR distribution of AGNs in our sample appears not to evolve with redshift. Despite 
X-ray AGNs and MS galaxies having different SFR distributions, the linear-mean SFR of 
AGNs derived from our distributions is remarkably consistent with that of MS galaxies, and 
thus with previous results derived from stacked Herschel data. This apparent contradiction 
is due to the linear-mean SFR being biased by bright outliers, and thus does not necessarily 
represent a true characterisation of the typical SFR of X-ray AGNs. 

Key words: galaxies: active—galaxies: evolution—galaxies: statistics 


1 INTRODUCTION 

Today’s most successful models of galaxy evolution predict that 
the energy released via accretion onto supermassive black holes 
(hereafter, BHs) has played an important role in dictating how to - 
day’s galaxies have grown and evolved (e.g.. lSchave et ahlboish . 
As such, understanding the connection between galaxy growth 
via star-formation and the growth of their resident BHs is one 
of the key challenges facing current extragalactic research. There 
are now numerous lines of empirical evidence in support of time- 
averaged/integrated BH growth correlating with star-formation in 


their host galaxies; for example, (a) the ti ght proportionality be - 
tween BH mass and galaxy bulge mass (e.g.. lGebhardt et al.l200(]h : 
(b) the similar cosmic histories of the volume-aver aged BH growth 
and star formation rat es (hereafter, SFR; e.g., ISilverman et al.l 
l2008l; lAird et al.l boish ; and, more directly, (c) the correlation 
between average BH g rowth and SFR among the star-forming 
galaxy population (e.g. , [Mullaney^taljj2()jjbh IChen et al.ir2013l ; 
iDelvecchio et al.1 1 20151 ; [Rodighiero et al. 20151) . However, it is 
still far from clear what physical processes (e.g., feedback pro¬ 
cesses/common fuel supply/common triggering mechanism) con¬ 
nect BH growth to star-formation to produce these average trends. 


* E-mail: j.mullaney@sheffield.ac.uk 


One of the primary means of making progress in this area 
has been to measure the SFRs and specific SFRs (i.e., SFR per 
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Figure 1. Host galaxy star-forming properties of our low-z (i.e., 0.5 < z < 1.5; not observed by ALMA) and high-z (i.e., z > 1.5) samples of AGNs (samples 
separated by the vertical dashed line). In all plots, grey circles indicate pre-ALMA (specific) star formation rates ([s]SFRs) from Herschel which are connected 
to their ALMA-measured (s)SFRs by dotted lines. (s)SFRs from ALMA are indicated by small white circles. Red and blue circles represent AGNs with Lx= 
JQ42-44 g-l Lx> 10"*^ ergs s“*, respectively, with lighter colours used for 3c upper limits. Top: SFR vs. redshift. Despite our ALMA observations 
probing SFRs up to a factor of fa 10 lower than Herschel, only fa 29 per cent of our ALMA-targeted AGNs are detected. Middle: sSFR vs. redshift. In this 
panel, the shaded region represents the average sSFR of main-sequence (MS) galaxies (SFRms) as described by Eqn. 9 of S 15 for the stellar mass range of our 
sample. Bottom: Rms vs. redshift. By definition, the horizontal line represents the average i?MS of MS galaxies. Shading indicates where i?MS< 0.5. Between 
34 and 55 per cent (dependent on upper limits) of AGNs in our combined (i.e., low-z-f high-z) sample lie within this shaded region, compared to Rf 15 per cent 
of MS galaxies. 


unit stellar mass, or sSFR) of galaxies hosting growing BFls (wit¬ 
nessed as active galactic nuclei, or AGN) and search for correla¬ 
tions or differences (vs. the non-AGN population) that may sig¬ 
nify a causal connection. The Herschel Space Observatory (here¬ 
after, Herschel) has played a major role in progressing this science 
by providing an obscuration-independent view of star-formation 
that is largely uncontaminated by emission from the AGN. Flow- 
ever, with even the deepest Herschel surveys detecting < 50 per 
cent of the AGN population, most studies have resorted to av er- 
aging (often via stacking analysis, but see lStanlev et al.ir2015h to 
characterise the (s)SFRs of the AGN population. These studies 
have typically reported that the average SFRs of AGNs trace that 
of star-forming “main - sequence” (hereafte r , MS) galaxies (e.g. , 
[ Mullaney et ^1 2012al; Santini et alj 20121 ; iHarrison et al ]M; 
[Rosario et al.ll2013l ; IStanlev et alJiTOlsI) . i.e.. the dominant popula¬ 
tion of star-forming galaxies whose SFRs are roughly proportional 
to their stellar mass (i.e., sSFR~constant), wi th a constant of pro¬ 
portionality that i ncreases with redshift (e.g.. iNoeske et al1l2007l ; 
iDaddi et al.ir2007h . Flowever, as averages can be biased by bright 
outliers, it is feasible that these finding s are being driven u pwards 
by a few bright sources (e.g.. Fig. 14 of iRosario et akll^OlSl) . Here, 
we test this by combining deep Herschel and ALMA observations 
to instead constrain the distribution of host galaxy SFRs of a sam¬ 
ple of X-ray selected AGNs and comparing it to that of MS galax¬ 


ies. We adopt Hq = 71 km s * Mpc = 0-73, = 0.27 and 

a Chabrier initial mass function (IMF). 


2 SAMPLE SELECTION 

To investigate any redshift evolution of the AGN (s)SFR distribu¬ 
tion, we use two samples of X-ray selected AGNs: a low-z sample 
spanning 0.5 ^ z < 1.5 and a high-z sample spanning 1.5 ^ z < 4 
(although the high-z sample is dominated by AGNs at 1.5 ^ z < 
2.7). The split at z = 1.5 is motivated by our ALMA target selec¬ 
tion criteria: for this, we only consider AGNs with redshifts >1.5 
since (a) the majority of z < 1.5 AGNs are detected with Herschel 
in the deepest fields and thus already have obscuration-independent 
SFR measures and (b) the negative k-correction at sub-mm wave¬ 
lengths would call for prohibitive ALMA integration times. 

The high-z sample were all selected from the 4 Ms Chandra 
Deep Field Sout h (hereafter, CDF-S) survey catalogue describ ed in 
IXue et all ( 1201 ih with updated redshifts from iHsu et alj | |2014|) ; we 
recalculate the rest-frame 2-10 keV luminosities (Lx) of the sources 
using these new redshifts. To ensure reliable AGN selection, we 
only consider those sources with Lx> 10^^ ergs s^* and reliable 
redshifts (spec-z, or phot-z with Az/(1 +z) < 0.1) that lie within 6' 
of the average aim point of the survey (the latter ensures reliable po- 
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Figure 2. The posterior probability distributions (PDs) for the parameters 
describing the assumed log-normal Rms distribution for AGN host galaxies: 
jd is the mode of the log-normal, while c is its la width (see Eq. 1). PDs for 
both our low-z and high-z samples are shown (see key). Contours of 20, 68 
and 95 per cent confidence are shown. The best-fit parameters of the com¬ 
bined (i.e., redshift-averaged) j?MS distribution of M S galaxies is indicated 
by the solid black circle tfrom ISchreiber et alJ201.5ll . The bottom and right¬ 
most plots indicate the relative probability of /j and a values; the location 
of the peak represent the most probable parameter values. When modelled 
as a log-normal, the Rms distribution of AGN host galaxies is significantly 
broader, and shifted significantly lower than that of MS galaxies. 


sitions for matching to ALMA counterparts). Our primary science 
goal of constraining the SFR distributions of AGN host galaxies in 
the context of the MS requires knowledge of the host galaxy s tellar 
masses (M*), which we derive following ISantini et al.l ll2012h . We 
refer to that study for a description of the relative uncertainties on 
M,, which is estimated to be 50 and 20 per cent (lo) for optically 
obscured (Type 2) and unobscured (Type 1) AGN, respectively. 
Since the majority (i.e., > 70 per cent) of the AGNs in our samples 
are optically obscured, this level of uncertainty has no significant 
impact on our conclusions. We restrict our sample to AGNs with 
M^:> 2 X 10^® Mq; below this threshold, it becomes prohibitive to 
reach low enough flux limits to probe to SFRs significantly below 
the mean SFR of MS galaxies (hereafter, SFRjvis) with ALMA. De¬ 
spite this Mt cut we still sample the vast majority of the luminous 
AGN population since the M* distri bution of Lx> 10^^ e rgs s^* 
AGNs peaks at 6 x 10^® Mq (e.g.. lMullanev et a02012al) . 

The above selection returned 49 AGNs (our high-z sample), 
with 20 and 29 having spec-z and phot-z, respectively. Of these 
49, 13 a re detected in the GOODS-Herschel 160 pm maps of the 
CDF-S dElbaz et aLllioTIh from which SFRs are be derived. Of 
the remaining 36 AGNs, 24 were observed by ALMA. Flowever, 
since making our original ALMA target list, a more sensitive Her- 
schel 160 pm map of the CDF- S has been generated by com- 
bining the PEP i Lutz et aTll201lh and GOODS-Herschel surveys 
dMagnelli et al.ll2013h and four of our 24 ALMA targets are now 
detected in that new map. For these four, we adopt the mean (s)SFR 
derived from the two facilities (see ij3). All other Herschel fluxes 
and 3 ct upper limits (including for the twelve Herschel-undetected 
AGNs not targeted by ALMA) are also taken from the combined 
PEP+GOODS-T/erjc/ie/ dataset. 


dra Deep Field North (from lAlexander et alj 

20031 and adopting 

the same redshifts and M* as iMullanev et al. 

2012a|) and South 


dXue etalJlioTl but using the updated redshifts and M*) sur- 


Table 1. Best-fit parameters for the log-normal Rms(= SFR/SFRms) dis¬ 
tributions (see Eqn. 1) of the samples of galaxies described in the main text. 


(1) 

(2) 

(3) 

Sample 

P 

G 

MS galaxies (Schreiber et al. 2015) 

-0.06“ 

0.31 ±0.02 

Low-z AGN sample 

-0.378+|>:0“ 

0.568«:“g 

High-z AGN sample 



Combined AGN sample 


0.560^°:«g 


Notes: Values given are the median of the posterior probability distribu¬ 
tions (PDs) and the 68 per cent confidence intervals. “This is slightly offset 
from exactly zero as Rms is the SFR relative to the linear mean SFR of MS 
galaxies, whereas p is the mode of the Rms distribution. 

veys with Herschel coverage by the PEP+GOODS surveys. We 
also restrict this low-z sample to Lx> 10"^^ ergs s^* and M,> 
2 X 10*** Mq to allow meaningful comparison with the high-z sam¬ 
ple. This returned a sample of 110 AGNs (i.e., our low-z sample), 
94 of which have spec-z. Sixty five of these 110 are detected in 
the Herschel 160 pm band, from which we derive (s)SFRs (see 
33; 3 ct flux upper limits were measured for the 45 Herschel non¬ 
detections. 


3 ALMA OBSERVATIONS AND DATA ANALYSIS 


All 24 of our ALMA targets were observed with ALMA Band- 
7 (i.e., observed-frame ~ 850 pm) during November, 2013, with 
a longest baseline of 1.3 km. To maximise observing efficiency, 
the ALMA-targeted sample was split into three groups according 
to the flux limit required to probe down to at least SERms at a 
given redshift. This corresponds to RMS flux limits of 200 pJy, 
125 pJy and 90 pJy for the three groups. ALMA continuum fluxes 
were measured using uv_f it of GILDAS v.aprl4c, adopting point 
source profiles for two unresolved sources and circular Gaussian 
profiles for the other five detected targets. 

Measured ALMA and Herschel fluxes and upper limits were 
converted to 8-1000 pm infrared luminosities (hereafter, Ljr) us¬ 
ing our adopted redshifts (see §2) and the ave r age in frared SEDs 
of MS galaxies described in IBethermin et al. i 2015h. which ar e 
constructed using the theoretical templates of IPraine & Lil (l2007h . 
However, we note that our conclusio ns do not change if we instead 
use either the lCharv & Elbaj (1200 ih SEDs or a starburst SED (i.e., 
Arp220). At the redshifts of our high-z sample, Band-7 probes the 
rest-frame 180-340 pm, close to the peak of the far-infrared emis¬ 
sion due to star-formation. Whil e these rest-frame w avelengths are 
also sens itive to dust mass (e .g., IScoville et al.ll201^ . based on the 
range of IPraine & Lil ( l2007h SED templates we estimate that the 
corresponding Ljr are accurate to within ±0.3 dex, which we fac¬ 
tor into our analyses. In a follow-up study we will employ full in¬ 
frared SED fitting incorporating all available Herschel and ALMA 
fluxes and upper limits to reduce the uncertainties associated with 
the adopted SED, but such detailed fitting is beyond the scope of 
this Letter. As a check, however, we note that the SERs derived 
from ALMA and Herschel data for the four AGNs that are de¬ 
tected with both are consistent to within this tolera nce. SERs are 
derived from Lir using Eqn. 4 from iKennicutJ (Il998h . but adopting 
a Chabrier IME. Einally, to explore the distributions of AGN host 
SFRs relative to SFRmSi we define Rjvjs = SFR/SFR^s, the rela- 
tive offset from the M S, where SFRms is computed using Eqn. 9 of 
ISchreiber et al.l ( l2015l ; hereafter, S15). 
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4 RESULTS 

4.1 Star-forming properties of X-ray AGNs 

Despite our ALMA observations probing to SFRs up to a factor of 
~ 10 below that achieved with Herschel (Fig. [T] top) only seven 
(i.e., 29 per cent) of the 24 ALMA-targeted AGNs in our high-z 

sample are detected at > 3o at 850 ;Um. The fractions of ALMA- 
undetected AGN are roughly the same for targets with spec-z and 
phot-z, suggesting that redshift uncertainties are not the primary 
cause of the non-detections. Despite the high fraction of non¬ 
detections, the 3 ct upper limits provided by the AhMA+Herschel 
data enable us to infer the level of consistency between the dis¬ 
tributions of Rms for AGN and MS galaxies (see ED, with the 
latter having been sh own not to vary in the and redshifts rang es 
considered here (e.e.. lRodighiero et'^l201 iLISargent et al.ll201 j) . 

To explore our AGN hosts’ star-forming properties in the con¬ 
text of the evolving MS, we plot their sSFRs and Rms values as a 
function of redshift (Fig.[T] middle and lower panels, respectively). 
We find that 54 to 88 (range due to upper limits) of the 159 AGNs 
(i.e., ~34 to ~55 per cent) in our combined (i.e., low-z-fhigh-z) 
sample have Rms< 0.5, with significant overlap between the frac¬ 
tions in our low-z (i.e., fii43 per cent to ~54 per cent) and high- 
z (i.e., ~14 per cent to ~59 per cent) samples. Comparing these 
fractions to the ;5il5 per cent of MS galaxies with Rms< 0-5 (from 
S15), reveals that the AGNs in our low-z sample, and possibly also 
our high-z sample, do not trace the same Rms distribution as MS 
galaxies, instead displaying a strong bias toward lower Rms values. 
Finally, we note that only ~5 per cent of AGNs in our combined 
sample reside in starbursts (i.e., with Rms> 4). 

4.2 Parameterlslng an X-ray AGN SFR distribution 

With the large fraction of AGNs with Rms< 0.5 in our combined 
and, in particular, low-z samples being inconsistent with the Rms 
distribution of MS galaxies, we now attempt to place constraints 
on the distribution of SFRs (relative to the MS; i.e., Rms) of AGN 
hosts. We place particular emphasis on quantifying the level of con¬ 
sistency/discrepancy between the AGN and MS Rms distributions. 

Our relatively small sample sizes, combined with the large 
fraction of non-detections prevents us from determining the AGN 
Rms distribution directly. Since a key goal here is to quantitatively 
compare the AGN and MS Rms distributions, we instead assume 
the same log-normal f orm for the AGN Rms distribution as found 
for MS galaxies fe.g.. lRodighiero et~^l201ll : ISargent et alJl2012L 
S15): 

(I) 

and infer its parameters (i.e., similar to lShao et al.ll20l5 who in¬ 
ferred the AGN Lir distribution). This is done purely to ease com¬ 
parison between the AGN and MS Rms distributions by allowing 
us to compare like-for-like parameters (i.e., the mode, and the 
variance, of the log-normal Rms distribution), and is not to be 
taken as a literal description of the true AGN Rms distribution^ 

We adopt a hierarchical Bayesian framework to determine 
the best-fit parameters (i.e., /j and o) for our assumed log-normal 
distributions, using Gibbs sampling and the Metropolis-Flastings 

* Investigating whether other forms better describe the Rms distribution of 
AGN hosts will be the focus of a later study incorporating a larger set of 
ALMA obseivations from Cycle 2 (PI: Alexander; awaiting completion). 



Figure 3. Rms distributions for our high-z and low-z samples of X-ray se¬ 
lected AGNs {Top) and MS galaxies (Bottom', from S15). Here, we show the 
log-normal distributions with best fitting parameters shown in Table 1 (solid 
and dotted curves; see key). The histograms in the top panel shows the rel¬ 
ative numbers of AGNs from our combined (i.e., low-z+high-z) sample in 
each Rms bin; the solid grey histogram represents those AGNs detected at 
> 3o with either Herschel or ALMA, whereas the empty histogram (with 
left-pointing arrows) also includes upper limits. The solid points in the top 
panel indicate the linear means of the log-normal distributions (equivalent 
to what would be obtained via, e.g., stacking analyses) and lie within lo of 
the linear mean Rms of MS galaxies (vertical dashed line). 

MCMC algorithm to rando mly sample their pos terior probability 
distributions (hereafter, PDs; lGelman et al1l2014h . The benefits of 
taking this approach are that (a) upper limits and uncertainties on 
Rms can he readily taken into account and (h) the resulting poste¬ 
rior PDs provide us with meaningful parameter uncertainties. We 
use weak prior PDs, noting that the centring of these priors (within 
reasonable limits) has no significant effect on our results. 

The posterior PDs on p and a for our two samples are pre¬ 
sented in Fig. [2 while the best-fit parameters (median of the PDs 
and 68 per cent confidence intervals) are given in Table[T] For com¬ 
parison, we also include the best-fit parameters of the log-normal 
Rms distribution for non-AGN MS galaxies from S15. As expected 
from the smaller size of our AGN sample and the high fractions of 
non-detections compared to the MS galaxy sample of S15, the un¬ 
certainties on the posterior parameter values for the assumed AGN 
log-normal Rms distribution are considerably larger than those for 
MS galaxies. Despite this, our analysis shows that the Rms distribu¬ 
tions of our low-z and high-z AGNs are both significantly broader 
and peak at significantly lower values (both at > 99.9 per cent con¬ 
fidence) than that of MS galaxies. Interestingly, our analyses show 
that, as with MS galaxies, there appears to be little evolution in the 
AGN Rms distribution, with the modes and variances of the log¬ 
normal distributions describing our low-z and high-z samples being 
consistent to within lo. In light of this, we infer the Rms distri¬ 
bution of our combined sample, which we find is roughly twice as 
broad as, and peaks ^ 0.4 dex below, that of MS galaxies (Table[Tll. 
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5 INTERPRETATION 

In the previous section we used our combined PdJsAA+Herschel 
data to demonstrate that, when modelled as a log-normal, the AGN 

distribution is significantly broader and peaks at significantly 
lower values than that of MS galaxies. This appears to be at-odds 
with recent findings based on mean-stacked Herschel data that the 
average star-forming propert ies of AGN hosts is c o nsistent with 
those of MS galaxies (e .g., iMullanev et al.l 12012al : ISantini et~n] 
l2012t[^sario et al.l2013h . Here, we place our results in the context 
of these studies to explore the root of these apparent discrepancies. 

When comparing to results derived from mean-stacked Her¬ 
schel data, it is important to note that mean-stacking provides a 
linear mean which will not correspond to the mode, /j, of a log¬ 
normal distribution. Instead, the linear mean will always be higher 
than the mode, with the discrepancy between the two increasing 
as function of both /j and <J. Therefore, while results from mean¬ 
stacking still hold when interpreted as the linear mean, depending 
on the underlying distribution this may not necessarily correspond 
to the mode. 

We compare our results against those from stacking by cal¬ 
culating the linear mean of our log-normal distributions, taking a 
Monte-Carlo approach to sample the /j and o PDs. This gives linear- 
mean AGNvalues (i.e., (/?ms)) of 0.99lQjg and 1-09^q25 
our low-z and high-z samples, respectively (Fig. O. These values 
are remarkably close to the linear mean of MS galaxies (i.e., 
{Rms)~ 1) and are broadly consistent with the linear means calcu¬ 
lated by mean-stacking Herschel 160 /rm maps at the positions of 
our AGN (i.e., (/?|y[s)= 0.81 ±0.12 and 0.86±0.15, respectively). 
We conclude that these linear-means are, indeed, influenced by the 
high tail of the broad distribution and may not necessarily give 
a reliable indication of the modal SFR of AGN hosts. 

Despite finding that the Rms distribution of AGN hosts is 
shifted toward lower values compared to MS galaxies, our results 
remain consistent with AGNs preferentially residing in galaxies 
with comparatively high (s)SFRs by z ~ 0 standards due to the 
strong redshift evolution of SFRms- Indeed, applying our analy¬ 
ses to sSFR (rather than Rms) gives distributions peaking at ~ 
0.2 Gyr^* and ^ 0.5 Gyr^* for our low-z and high-z samples, re¬ 
spectively. To put this in context, (sSFRms) ^0.1 Gyr^* at z ~ 0, 
thus local galaxies with sSFRs of 0.2 Gyr^* and 0.5 Gyr^* would 
be classed as MS and starbursting galaxies, respectively. 

Our result compare favourably to those derived from AGN sur¬ 
veys conducted at other wavele ngths. For examp l e, usin g SFRs de¬ 
rived from optical SED fitting, iBongiomo et alJ l l2012h reported a 
broad sSFR distribution for X-ray selected AGNs that peaks be¬ 
low that of the MS at red shifts similar to tho se explored here (i.e., 
0.3 < z < 2.5). Similarly, I Azadi et al.l ( I20l3) showed that the Rms 
distribution of X-ray selected AGNs (with a similar M* selection 
as here) peaks at ~ 0.1 and is similar to the Rms distribution of 
M*-matched galaxies (i.e., not just star-forming galaxies). As such, 
these studies and the results presented here support the view that 
X-ray selected AGN hosts at moderate to high redshifts span the 
full range of relative sSFRs of 2 x 1 0*° M© galaxies (e.g., 
iBrusa et^l2009l : iGeorgakakis et al.ll2014h . However, with recent 
results suggesting that X-ray a bsorbed AGN may have higher SFRs 
than unabsorbed AGN (e.g.. I Juneau et alJ 12013l : lOel Moro etakl 
|20H . it is feasible that alternative AGN selections may bring the 
AGN Rms distribution closer to that of MS galaxies. 

We thank the anonymous referee. DMA, ADM, CMH acknowledge 
STFC grant ST/1001573/1. This paper makes use of ALMA data: 
ADS/1AO.ALMA#2012.1.00869.S. 
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